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The general equations of motion for ocean dynamics are presented and the waves supported 
by the (inviscid, unforced) linearized system with respect to a state of rest are derived. The 
linearized dynamics sustains one zero frequency mode (called buoyancy mode) in which salinity and 
temperature rearrange in such a way that seawater density does not change. Five nonzero frequency 
modes (two acoustic modes, two inertia-gravity or Poincare modes, and one planetary or Rossby 
mode) are also sustained by the linearized dynamics, which satisfy an asymptotic general dispersion 
relation. The most usual approximations made in physical oceanography (namely incompressibility, 
Boussinesq, hydrostatic, and quasigeostrophic) are also consider, and their implications in the 
reduction of degrees of freedom (number of independent dynamical fields or prognostic equations) 
of, and compatible waves with, the linearized governing equations are particularly discussed and 
emphasized. 

PACS numbers: 43.30.Bp, 43.30.Cq, 43.30.Ft 



I. INTRODUCTION 

The goal of this educational work is to show how the 
various types of linear waves in the ocean (acoustic, 
inertia-gravity or Poincare, and planetary or Rossby 
waves) can be obtained from a general dispersion relation 
in an approximate (asymptotic) sense. Knowledge of the 
theory of partial differential equations, and basic classical 
and fluid mechanics are only needed for the reader to 
understand the material presented here, which could be 
taught as a special topic in a course of fluid mechanics 
for physicists. 

The exposition starts by presenting the general 
equations of motion for ocean dynamics in Sec. [H] This 
presentation is not intended to be rigorous, but rather 
conceptual. Accordingly, the equations of motion are 
simplified as much as possible for didactic purposes. The 
general dispersion relation for the waves supported by 
the (inviscid, unforced) linearized dynamics with respect 
to a quiescent state is then derived in Sec. IIIII This is 
done by performing a separation of variables between the 
vertical direction, on one side, and the horizontal position 
and time, on the other side. Sec. IfVI discusses the 
implications of the most common approximations made 
in oceanography (namely incompressibility, Boussinesq, 
hydrostatic, and quasigeostrophic) in the reduction of 
degrees of freedom (number of independent dynamical 
fields or prognostic equations) of, and compatible waves 
with, the linearized governing equations. Particular 
emphasis is made on this important issue, whi ch i s 
vaguely covered in standard textbooks (e.g. Refs. M llol) 
Some problems have been interspersed within the text to 
help the reader to assimilate the material presented. The 
solutions to some of these problems are outlined in App. 
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II. GENERAL EQUATIONS OF MOTION 

Let x := (x,y) be the horizontal position, i.e. 
tangential to the Earth's surface, with x and y its 
eastward and northward components, respectively; let z 
be the upward coordinate; and let t be time. Unless 
otherwise stated all variables are functions of (x,z,i) in 
this paper. 

The thermodynamic state of the ocean is determined by 
three variables, typically S (salinity), T (temperature), 
and p (pressure, referred to one atmosphere). Seawater 
density, p, is a function of these three variables, i.e. p = 
p{S,T,p), known as the state equation of seawater. In 
particular, 

p- 1 Bp = a s BS-a T BT + a p Bp. (1) 

Here, D := dt + u • V + wd z is the substantial or 
material derivative, where u and w are the horizontal and 
vertical components of the velocity field, respectively, and 
V denotes the horizontal gradient; as ■= p^ 1 {dsp) Tp 
and ut '■= p i^Tp)g v are the haline contraction and 
thermal expansion coefficients, respectively; and a p := 
p^ 1 (d p p) T s = o>tT + p^CjT 2 , where T is the adiabatic 
gradient and c s is the speed of sound, which characterize 
the compressibility of seawater. 

The physical state of the ocean is determined at every 
instant by the above three variables (S,T,p) and the 
three components of the velocity field (u, w), i.e. six 
independent scalar variables. The evolution of these 
variables are controlled by 

BS = F S , (2a) 

BT = TBp + F T , (2b) 

Bp = -a.- 1 (V • u + d z w + a T F T - a s F s ) , (2c) 

Du=-p- 1 Vp + F u , (2d) 

Bw = -p- 1 d zP -g + F w . (2e) 

In Newton's horizontal equation (0i), the term F u 
represents the acceleration due to the horizontal 
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components of the Coriolis and frictional forces. In 
Newton's vertical equation the term F w represents 
the acceleration due to the vertical component of the 
Coriolis and frictional forces, and g is the (constant) 
acceleration due to gravity. The term F$ in the salinity 
equation J2t) represents diffusive processes. The term 
Ft in the thermal energy equation {2b) j which follows 
from the first principle of thermodynamics, represents 
the exchange of heat by conduction and radiation, as 
well as heating by change of phase, chemical reactions or 
viscous dissipation. The pressure or continuity equation 
follows from {TJ|. 

Problem 1 Investigate why (01, e) do not include the 
centrifugal force which would also be needed to describe 
the dynamics in a noninertial reference frame such as one 
attached to the Earth. 

Since adiabatic compression does not have important 
dynamical effects, in physical oceanography it is 
commonly neglected. This is accomplished upon 
introducing the potential temperature 9, which satisfies 
a e B9 = a T (DT - TDp), so that p = p(S, 9,p) and QJ is 
consistently replaced by 

p- 1 Bp = a s BS-a e m + p- 1 c; 2 Bp, (3) 

where here it must be understood that as = {dsp)e.p 
and a-r = (dgp)g tP . Equations (03, c) then are replaced, 
respectively, by 

D9 = F g , (4a) 
T>p=-pc 2 s {V -u + dzW + aeFe - a s F s ) . (4b) 

As our interest is in the waves sustained by the 
linearized dynamics, we do not need to consider cither 
diffusive processes or allow the motion to depart from 
isentropic. Hence, we will set Fs = = Fg so that 
equations an d can be substituted, respectively, 
by 

DC = w, (5a) 
Dp = -pc 2 s (V • u + d z w) , (5b) 

where C is the vertical displacement of an isopycnal which 
is defined such that p — p T (z — £). 

Problem 2 Show that equations {2Jl) and certainly 
lead to {30 when F s = = Fg . 

We will also neglect frictional effects, so that equations 
{2Ji,e) are seen to be nothing but Euler equations of 
(ideal) fluid mechanics with the addition of the Coriolis 
force. The latter will be further considered as due solely 
to the vertical component of the Earth rotation. Thus, 
the following simplified form of equations {2JI,e) will be 
considered: 

Du = — p^ 1 Vp — fz x u, (6a) 
Dw = -p~ 1 d z p - g. (6b) 



Here, z is the upward unit vector and / := 20 sin 
where is the (assumed constant) spinning rate of the 
Earth around its axis and •d is the geographical latitude, 
is the Coriolis parameter. For simplicity, we will avoid 
working in full spherical geometry. Instead, we will 
consider / = /o + fly, where /q := 2fisini?o and '■= 
2fii? _1 cosi?o with -do a fixed latitude and R the mean 
radius of the planet, and V = (d x ,d y ), which is known 
as the /3-plane approximation. It should remain clear, 
however, that a consistent /3-plane approximation must 
include some geometric (non-Cartesian) termsi Neither 
these terms nor those of the Coriolis force due to the 
horizontal component of the Earth's rotation contribute 
to add waves to the linearized equations of motion. Their 
neglection is thus well justified for the purposes of this 
paper. 



III. WAVES OF THE LINEARIZED DYNAMICS 

Consider a state of rest (u = 0, w = 0) characterized 
by dp r /dz = —p T g, where p r {z) and p T (z) are reference 
profiles of pressure and density, respectively. Equations 
Q and ©, linearized with respect to that state, can be 
written as 



dt(' = w', (7a) 

d t p' = -p t cl (V • u' + d-w') , (7b) 

d t u' = -p- 1 V P ' - fixu', (7c) 

d t w' = -p^O+p' - N 2 ('. (7d) 



Here, primed quantities denote perturbations with 
respect to the state of rest; df := d z ± gc~ 2 ; N 2 (z) := 
—g(p~ 1 dp 1 ./dz + gc~ 2 ) is the square of the reference 
Brunt- Vaisala, frequency; and c s is assumed constant. In 
addition to the above equations, it is clear that 

d t S' = -wdSr/dz, (8a) 
d t 9' = -wd6 T /dz, (8b) 

where S r (z) and 9 r (z) are reference salinity and potential 
density profiles, respectively. 

Problem 3 Work out the linearization of the equations 
of motion. 

A. Zero Frequency Mode 

The linearized dynamics supports a solution with 
vanishing frequency (dt = 0) such that 

C' EE 0, p' EE 0, U' = 0, w' EE 0, (9) 

as it follows from {7J|, but with 

sVo, e'^o, (io) 



3 



as can be inferred from (JSJ). Namely for this solution 
the salinity and temperature fields vary without changing 
the density of the fluid. More precisely, one has, on one 
hand, p' = p 1 .(g^ 1 N 2 + gc~ 2 )£ , and, on the other, p' — 
p T (asS' — ag8' + a p p'), where the a's are evaluated at 
the reference state. By virtue of © then it follows that 



asS' — agO' = 0. 



(11) 



This so-called buoyancy mode describes small scale 
processes in the ocean such as double diffusion. 



B. Nonzero Frequency Modes 

Upon eliminating w between (JJJi) and (01), and 
proposing a separation of variables between z, on one 
side, and (x, t), on the other side, for the horizontal 
velocity and pressure fields in the form 

u' = u c (x,t)9 z -F(z), p' = Pr (z)p c (x,t)9;F(z), (12) 

it follows that 

c: 2 d;F8 tP c + (d~FV ■ u c + d- t £) = 0, (13a) 
d t u c + fz x u c + Vp c = 0, (13b) 
duC' + P^dt { Pr d-F) p c + N 2 £ = 0. (13c) 

Now, assuming a common temporal dependence of the 
form e~ luJt , from Q13b) one obtains 



c' = 



P c d+ ( Pl d-F) 

p T (N 2 -uj 2 ) ' 



(14) 



Then, upon substituting (jl4|) in l|13h ) it follows that 



d7F 



-d: 



d+ ( Pl d-F) 



p T {N 2 - uj 2 ) 



V ■ u c 

d t p c 



= cr 2 , (15) 



where c is a constant known as the separation constant. 
Clearly, c s must be chosen as a constant in order for 
the separation of variables to be possible. From 115|) it 
follows, on one hand, that 

dt ( Pl d-F) + Pr (N 2 - u; 2 ) (c~ 2 - c- 2 ) F = 0, (16) 

and taking into account (113b ). it follows, on the other 
hand, that 



d t p c + c 2 V -u c = 0, 
d t u c + fz x u c + Vp c = 0. 



(17a) 
(17b) 



Equation 1)16(1 governs the vertical structure of the 
perturbations, whereas system H17|) controls the evolution 
of their horizontal structure. 

Problem 4 Show that the alternate separation of 
variables which uses F(z) instead of d~F(z) leads to 
a vertical structure equation with a singularity where 
uj 2 = N 2 . 



Equation l|16|) can be presented in different forms 
according to the approximation performed. Under 
the incompressibility approximation, which consists of 
making the replacement d t p' + gw' i— > in the continuity 
equation J3d), equation l|16|) takes the form 



dt ( Pl d z F) + Pl c~ 2 (N 2 ~u 2 )F = 0. 



(18) 



This approximation corresponds formally to taking the 



limit 



0. The hydrostatic approximation, in turn, 



consists of making the replacement d t w' i— ► in Newton's 
vertical equation (|7Ji). This way, without the need 
of assuming any particular temporal dependence, it 
follows that (' — —p c d+ (p r d~ F) /(p r N 2 ). Consequently, 
equation p6|l reduces to 



d+(p T d-F)+ Pl (c- 



c~ 2 ) N 2 F = 0. (19) 



This approximation is valid for uj 2 <C N 2 , i.e. periods 
exceeding the local buoyancy period which typically is of 
about 1 h. (Of course, this approximation implies that 
of incompressibility as it filters out the acoustic modes 
whose frequencies are much higher than the Brunt- 
Vaisala, frequency.) Another common approximation is 
the Boussinesq approximation, which consists of making 
the replacements p r h- > p — const, and >— > d z in J7|. 
Under this approximation, equation takes the simpler 
form 



d 2 F/dz 2 + (c- 2 - c- 2 ) (N 2 -uo 2 )F = 0. 



(20) 



Problem 5 Show that the Boussinesq approximation is 
very good for the ocean but not so for the atmosphere. 
Hint: This approximation requires c 2 3> gH, where H is 
a typical vertical length scale. 



1. Horizontal Structure 

To describe these waves is convenient to introduce a 
potential </?(x, t) such tha^SiS 



P C = -C 2 (d t y + fd X ) if, 

u c = (c 2 d xy + fdt) V , 
v c = (d tt - c 2 d xx ) ip, 



(21a) 
(21b) 
(21c) 



which allows one to reduce system <|17[) to a single 
equation in one variable: 

:= {d t [d tt + f 2 (v) - c 2 V 2 ] - /3c 2 d x } ip = 0. (22) 

The linear differential operator £ contains a variable 
coefficient and, hence, a solution to (|22|l must be of the 
form 

ip = $(y)e i(fca; - wt) (23) 
with $(y) satisfying 

d 2 $/dy 2 + l 2 (y)<P = (24) 
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where 



f(y):=-k*-^ + ^-pyl. (25) 

Now, if l 2 (y) is positive and sufficiently large, then Q(y) 
oscillates like 



(26) 



This is known as the WKB approximation (cf. e.g. Ref. 
4), where l(y) defines a local meridional wavenumber in 
the approximate (asymptotic) dispersion relation 



u; 2 -(f 2 +c 2 k 2 )-f3-=0, 



(27) 



where k := (k, I) is the horizontal wavenumber. 

System ljl 71) also supports a type of nondispersive 
waves called Kelvin waves. These waves have v c = 
and thus are seen to satisfy 



d t p c + c 2 x u c = 0, 

d t u c + d xP c = o, 
fu c + d yP c = o. 



(28a) 
(28b) 
(28c) 



Clearly, these waves propagate as nondispersive waves in 
the zonal (east-west) direction — as if it were / = — 
and are in geostrophic balance between the Coriolis and 
pressure gradient forces in the meridional (south-north) 
direction. From l|28b . b) it follows that 



p c = A{y)K(x - ct) = cu c , 



(29) 



where K (•) is an arbitrary function. By virtue of (|28b ~) 
then it follows dA/dy + fA/c = 0, whose solution is 

A(y) oc e-T d y f^/ c = e -(/^+5^ 2 )/ Cj 

which requires, except there where /q = (i.e. the 
equator), the presence of a zonal coast to be physically 
meaningful. 

Problem 6 Consider the Kelvin waves in the so-called 
/ plane, i.e. with (3 = 0. 



2. Vertical Structure 

Under the Boussinesq approximation the five fields 
of system Q remain independent, thereby removing no 
wave solutions. We can thus safely consider the vertical 
structure equation l|2U|) . which we rewrite in the form 



where 



d 2 F/dF 2 +m 2 {z)F = 



m 2 (z) := [N 2 {z)~lo 2 ] (c 



(30) 



(31) 



Equation l|30|) can be understood in two different senses. 
Within the realm of the WKB approximation, l|31|) 
defines a local vertical wavenumber, and a solution to 
(|30|l oscillates like 



F(z) 



^ii J dz m(z) 



(32) 



The other sense is that of vertical normal modes, in which 
(|30|l is solved in the whole water column with boundary 
conditions 

F(-H) = 0, (33a) 
g F(0) = c 2 dF(0)/dz. (33b) 

Condition l|33b l comes from imposing w 1 — at z = — H 
where H, which must be a constant, is the depth of the 
fluid in the reference state. Condition \6'6b ) comes from 
the fact that p' = at z = 0, which means that the 
surface is isopycnic (i.e. the density does not change 
on the surface). This way one is left with a classic 
Sturm-Liouville problem. Making the incompressibility 
approximation and assuming a uniform stratification in 
the reference state, namely N — N = const., it follows 
that 



lo 2 = N 2 



rag tan mH. 



(34) 



(Notice that to obtain m is necessary to fix a value of 
uj.) In the hydrostatic limit, lo 2 -C ^V 2 , it follows that the 
vertical normal modes result from 

tan mH = s/(mH), (35) 

where s := N 2 H/g, which is a measure of the 
stratification, is such that < s < oo by static stability. 
In the ocean s is typically very small, so for s <g; 1 from 
|jS5J it follows that 



N/yfg~H if % = 0, 
iir/H if i = 1,2, 



(36) 



The first mode is called the external or barotropic mode; 
the rest of the modes are termed the internal or baroclinic 
modes, which are well separated from the latter in what 
length scale respects. More precisely, the Rossby radii 
of deformation are defined by Ri := N/(rrij |/o|); for 
the barotropic mode Ro = \/gH/\fo\ whereas for the 
baroclinic modes Ri = NH/(iir \fo\) = y/sRo/(i7r) <C 
Ro. Finally, the rigid lid approximation consists of 
making w' — at z = 0, which formally corresponds to 
take the limit g — > oo in (|33b ) . This approximation filters 
out the barotropic mode since it leads to tan mH = 0. 

Problem 7 Demonstrate that p' = gC? at z = 0. 



3. General Dispersion Relation 

Upon eliminating c between (I25|) and (|31|l it follows 
that 



k 2 + (3k /ui m 2 
uj 2 - P = N 2 -uj 2 



(37) 
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which is a fifth-order polynomial in ui that constitutes 
the general dispersion relation for linear ocean waves in 
an asymptotic WKB sense. This is the main result of 
this paper. Approximate roots of l|37l) are: 

acoustic : J 2 = (k 2 + m 2 )c 2 , (38) 

which holds for lj 2 ^> iV 2 (i.e. very high frequencies); 



Poincare : ui = 



, k 2 7V 2 + m 2 / 2 



(39) 



which follows upon taking the limit c s 2 —> and is valid 
for frequencies in the range f 2 <ui 2 < N 2 ; and 



Rossby : u> 



8k 



k 2 + {m 2 /N 2 )f 2 ' 



(40) 



which also follows in the limit c s 2 — > but is valid for 
lo 2 <C f 2 (i.e. very low frequencies). 

Problem 8 Demonstrate that the classical dispersion 
relations for Poincare waves, lo 2 = f 2 + c 2 k 2 , and surface 
gravity waves, lo 2 — g |k| tanh |k| H, are limiting cases of 

CD. 



IV. DISCUSSION 




2 AW 

5 independent fields <-> 5 waves : < 2 PW 

1 RW 



incompressibility J, 



d t C = ■ 

d t u' = 



d z w' = -V • u' 

A P ' = -V • fz x u' - <9 Z (7V 2 

' 2 PW 
1 RW 



3 independent fields <-> 3 waves 



iV 

quasigeostrophy j c? z u = — — z x V£ 



fin 7 



/o 



1 independent field «-» 1 wave : 1 RW 



TABLE I: Reduction of independent fields (and, hence, 
prognostic equations) by the incompressibility and quasi- 
geostrophic approximations. Here, AW, PW, and RW 
stand for acosutic waves, Poincare waves, and Rossby wave, 
respectively; A := V 2 + d zz is the three-dimensional 
Laplacian; and q := f + V 2 p'//o + d z (foN~ 2 d z p') is the 
so-called quasigeostrophic potential vorticity. 



The inviscid, unforced linearized equations of motion 
(|7|l have five prognostic equations for five independent 
dynamical fields. As a consequence, five is the 
number of waves sustained by (J7J| which satisfy the 
general dispersion relation (|37|l in an asymptotic WKB 
sense. In proper limits, two acoustic waves (AW), 
two Poincare waves (PW), and one Rossby wave (RW) 
can be identified. The fact that the number of waves 
supported by the linearized dynamics equals the number 
of independent dynamical fields or prognostic equations, 
i.e. the degrees of freedom of (JZJ), means that the waves 
constitute a complete set of solutions of the linearized 
dynamics (cf. Table |]J. 

The number of possible eigensolutions can be reduced 
is approximations that eliminate some of the prognostic 
equations, or independent dynamical fields, of the system 
are performed. The Boussinesq approximation, which 
is very appropriate for the ocean, does not eliminate 
prognostic equations and has the virtue of reducing 
the mathematical complexity of the governing equations 
considerably. The incompressibility approximation, in 
turn, removes two degrees of freedom: the vertical 
velocity is diagnosed by the horizontal velocity, 

cW = -V • u', (41) 

and the latter along with the density diagnose the 
pressure field through the three-dimensional Poisson 
equation 

(V 2 + d zz )p' = -V • fz x u' - d z (N 2 ('). (42) 



As a consequence, the two AW are filtered out and 
one is left with the two PW and the RW. With 
these two approximations the Euler equations Q 
reduces to the so called Euler- Boussinesq equations. 
Performing in addition the hydrostatic approximation, 
which corresponds to neglecting dtw' in Newton's vertical 
equation, does not amount to a reduction of degrees of 
freedom because the vertical velocity is already diagnosed 
by the horizontal velocity. In this case, the density field 
diagnoses the pressure field through 

d zP ' = -N 2 ('. (43) 

With this approximation (which implies that of incom- 
pressibility) and the Boussinesq approximation, system 
reduces to what is known in geophysical fluid 
dynamics as the primitive equations. Finally, one 
approximation that eliminates independent fields is the 
quasigeostrophic approximation, which is often used to 
study low frequency motions in the ocean, and the Earth 
and planetary atmospheres. In this approximation the 
density diagnoses the horizontal velocity through the 
"thermal wind balance," 

_/V 2 

d z W = — z x VC', (44) 
Jo 

thereby removing two degrees of freedom and leaving 
only one RW. 
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APPENDIX A: SOLUTIONS TO SOME OF THE 
PROBLEMS 

Problem 1 To describe the dynamics in a noninertial 
reference frame such as one tied to the rotating Earth, 
two forces must be included: the Coriolis and centrifugal 
forces. However, Laplacei showed that if the upward 
coordinate z is chosen not to lie in the direction of the 
gravitational attraction, but rather to be slightly tilted 
toward the nearest pole, the centrifugal and gravitational 
forces can be made to balance one another in a horizontal 
plane (cf. also Ref. |8|). With this choice the Coriolis force 
is the only one needed to describe the dynamics. Notice 
that the absence of the centrifugal force in a system fixed 
to the Earth is what actually makes rotation effects real: 
they cannot be removed by a change of coordinates. 

Problem 2 In the absence of diffusive processes, the 
isopycnal z = ^ is a material surface, i.e. [u + z w — 
(u c + zw c )} ■ [VC - 2(1 - dgQ] = 0. Here, u c + zw c 
denotes the velocity of some point on the surface [the 
velocity of a surface is not defined and it only makes 
sense to speak of the velocity in a given direction, e.g. 
the normal direction, in whose case it is z (1 — d z Q — V£]. 
From the trivial relation z — £ = then it follows that 
(u<j + zii)()-[V(-z (1— d z Q] = —d t ( and, hence, DC = w 
at z — C- 

Problem 3 To perform the linearization of the equa- 
tions of motion, we write 

(u,w,C) = (u>',0 + ••• , 

(p,p) = (p t ,pr) + {p',p') + ■■■ , (Al) 

0:1 a a 2 

where a is an infinitesimal amplitude. The 0(a) 



continuity equation J7|d) readily follows upon noticing 
that, up to O(a), c~ 2 Dp — c~ 2 (d t p' — p T gw'). Up to O(a), 
Dp — c~ 2 Dp = d t p' — g~ 1 p 1 .N 2 w' — c~ 2 d t p' and DC — w = 
dt(' — w'. Then from the relationships Dp — c~ 2 T)p = 
and DC - w = it follows that p' = c~ 2 p' + g~ 1 p ] .N 2 ('. 
Bearing in mind the latter relation and the fact that 
dp r /dz = —gp r , the 0(a) vertical Newton's equation (JZJi) 
then follows. 

Problem 4 For the ocean c s ~ 1500 m s^ 1 3> \fg~H ~ 

200 m s~ x ; by contrast, for the atmosphere c s ~ 350 m 
s _1 ~ \/ r gH with H ~ 12 km, which is the typical height 
of the troposphere. 

Problem 7 At the surface z — r] it is w = d t r\ + u • Vrj 
and p — (here, p is a kinematic pressure, i.e. divided by 
a constant reference density p). Writing r\ = r( + 0(a 2 ) 
and Taylor expanding about z = it follows, on one 
hand, 

w' + rj'd z w' + 0{a 3 ) = d t r{ + u' • Vr/ + 0(a 3 ) (A2) 

at z = 0, and, on the other hand, 

Pl + (p r + dp r /dz)r/ + V 'd zP ' + 0(a 3 ) = (A3) 

at z = 0. From i|A21 it follows, to the lowest order, w' = 
dtf)' at z = 0. Since w' = dt(' , for a wave (i.e. dt ^ 0) 
then it follows that r/' — C' at z = 0. Taking into account 
the latter and choosing p = /9 r (0), from l|A3|l it follows, 
to the lowest order, p' = gr/' = gQ' at z = since p r = 
and dp r /dz = —gp r /p = —g at z = 0. 

Problem 8 The classical dispersion relation for 
Poincare waves corresponds to the hydrostatic limit, 
which requires m 2 ^S> k 2 (i.e. that the vertical length 
scales be shorter than the horizontal length scales). 
Under this conditions, uj 2 — f 2 + k 2 N 2 /m 2 = f 2 + c 2 k 2 . 
To obtain the dispersion relation for surface gravity 
waves one needs to take into account boundary conditions 
jjsSjJl : making TV 2 = = f 2 it follows, on one hand, that 
m 2 = — k 2 , and, on other, that uj 2 = —mgt&YimH. The 
dispersion relation uj 2 = g |k| tanh |k| H then readily 
follows. 
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